Preservation of quantum correlation between separated nitrogen-vacancy centers 

embedded in photonic crystal cavities 
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We investigate the non-Markovian dynamics of quantum correlation between two initially entan- 
gled nitrogen- vacancy centers (NVC) embedded in photonic crystal cavities (PCC). We find that 
a finite quantum correlation is preserved even asymptotically when the transition frequency of the 
NVC is within the band gap of the PCC, which is quantitatively different from the result of ap- 
proaching zero under the Born-Markovian approximation. In addition, once the transition frequency 
of NVC is far beyond the bad gap of the PCC, the quantum correlation initially prepared in NVC 
will be fully transferred to the reservoirs in the long-time limit. Our result reveals that the inter- 
play between the non-Markovian effect of the structured reservoirs and the existence of emitter-held 
bound state plays an essential role in such quantum correlation preservation. This feature may open 
new perspectives for devising active decoherence-immune solid-state optical devices. 
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I. INTRODUCTION 

Recently, decoherence control and non-Markovian dy- 
namics due to strong backaction from the environment 
have attracted much attention in practical implementa- 
tion of nanoscale solid-state quantum information pro- 
cessing (QIP) both theoretically [1-4] and experimen- 
tally [5-9]. It is widely believed that the unavoidable 
interaction with Markovian environment results in de- 
coherence effect and loss of quantum correlation, so the 
quantum correlation preservation (QCP) in quantum sys- 
tems has become a critical challenge in quantum compu- 
tation technologies. Previous studies have shown that 
long-time quantum entanglement protection can be re- 
alized by engineering the structured environment such 
as photonic band gap (PBG) materials, which has peri- 
odic dielectric structures exhibiting a range of frequencies 
with electromagnetic wave propagation forbidden [10]. 
This novel property mostly arises from the specific struc- 
ture of PBG environment [11], which leads to strong 
emitter-field correlation and formation of emitter-field 
bound states (EFBS) [12] with the spontaneous emission 
being greatly inhibited [13]. Experimentally, the non- 
Markovian effect of the emitter in the PBG reservoir has 
been directly observed in several kinds of nanoscale cav- 
ities [7-9]. 

In this Letter, we focus on the issue of decoherence 
suppression in solid-state systems consisting of diamond 
nitrogen-vacancy centers (NVCs) [14] and planar pho- 
tonic crystal cavities (PCC) [15]. The NVC is an attrac- 
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FIG. 1: (Color online) The composite NV-PCC system con- 
sists of a planar PCC and two identical NVCs in diamond 
nanocrystals. The inset shows the level structure of a NVC, 
where the electronic ground state | 3 ^4) is a spin triplet state, 
and D eg = 7 e Bo is the level splitting induced by an external 
magnetic field Bo with 7 e the electron gyromagnetic ratio. 
The red arrow denotes the coupling between NVC and PCC. 
r is the spontaneous decay rate of the excited state I-A2). 



tive spin qubit since it exhibits a unique combination 
of robust room-temperature spin coherence [16] and ef- 
ficient optical addressability, controllability, and readout 
[17-20]. Furthermore, the planar PCC with high Q fac- 
tor (i.e., Q > 1 x 10 6 even up to 2 x 10 7 in the optimized 
structure [21, 22]) and mode volumes comparable to a cu- 
bicoptical wavelength [23] can strongly confine photons 
in a tiny space of optical-wavelength dimension within 
a PBG structure. Numerous experiments have succes- 
sively demonstrated the strong coupling between NVCs 
and the modes of silicon nitride PCC [24], gallium phos- 
phide PCC [21, 23, 25], and PCC in monocrystalline di- 
amond [26, 27], respectively. These advances imply that 
the conventional Born-Markovian description to the NVC 
decoherence induced by the PPC is not applicable any- 
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more. 

In the present work we investigate the non-Markovian 
dynamics of quantum correlation between two initially 
entangled NVCs, each of which is embedded in a 
nanocavity in the planar PCC coupled to the radiation 
fields initially in vacuum states. Characterizing the quan- 
tum correlation by both quantum discord (QD) [28] and 
entanglement of formation (EoF) [29], we find that, if 
the transition frequency of the NVC is fully within the 
band gap of the PCC, a finite quantum correlation is pre- 
served due to the interplay between the non-Markovian 
effect of the structured reservoirs and the existence of 
EFBS. Otherwise, no QCP can be observed, and there 
only exists the quantum correlation between the two in- 
dependent reservoirs in the long-time limit. The essential 
condition for realizing this QCP is explicitly given. These 
results would be useful for experimental exploration of 
non-Markovian features in spin-based quantum system. 



II. THE MODEL AND HAMILTONIAN 

Now we focus on two noninteracting NVCs (NV1 and 
NV2) coupled, respectively, to two uncorrelated vacuum 
reservoirs R\ and i?2, namely, two nanocavities in a pla- 
nar PCC, as shown in Fig. 1. Each NVC is nega- 
tively charged with two unpaired electrons located near 
the vacancy, usually treated as electron spin-1. In our 
scheme, the PCC modes with a + polarization are cou- 
pled to the transition from the NVCs ground state sub- 
level |— ) = | 3 A,m s = — f) to one of the excited states 
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+l) + \E+,m B = -l))/V2[20]. 



Due to the dynamical independence between the two lo- 
cal subsystems, we can first solve the individual subsys- 
tem, and then apply the result obtained to the double 
case. Each subsystem (i.e., a NVC plus a radiation field 
propagating in the PCC) is governed by [30] 



w k a\a k 



+ 22 k (gkV+a k + h.c.), (1) 



where a± and loq are the inversion operator and transi- 
tion frequency of the NVC, and a\ (a k ) are the creation 
(annihilation) operators of the fc-th mode of the reservoir. 
The coupling strength between the NVC and the reser- 
voir is denoted by g k = uJad ■ e& / V2eoWfcTo [31], where 
e k and Vo are the unit polarization vector and the nor- 
malization volume of the radiation field, d is the dipolc 
moment of the NVC, and eo is the free space permittivity. 
Here the specific periodic structure of the PCC causes a 
band-gap dispersion relation to the field [32]. 



uj k = uj c + A(k - k y 



(2) 



where A ~ w c /fcg with fco ~ uj c /c being a specific wave 
vector with respect to the point-group symmetry of the 
PCC, and uj c is the dip of the band frequency. Note that 
Eq. (2) describes well our 2D PPC structure, where the 
2D square lattice is formed by cylinders (see Fig. 1). 



Here the two orthogonal directions are equivalent and 
the dispersion relation is the same for both. 
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FIG. 2: (Color online) Time evolution of the quantum corre- 
lations between NV1 and NV2 calculated by QD (solid-line) 
and EoF (dashed-line), where (a) a = l/y/2, and ujq = ui c /W 
(b) a = 0.2, and wo = w c /10; (c) a = 0.3, and ljq = w c /W 
(d) q = l/y/2, and ujq = Wuj c ; (e) a = 0.2, and ujq = 10co c 
(f) a = 0.3, and luq = 10lu c . The parameter rj — 0.2 is used. 



III. THE DYNAMICS OF QUANTUM 
CORRELATIONS 



We firstly study the time evolution of a single NVC 
prepared initially in the excited state |+) influenced by 
the reservoir. The state of the system evolves as 

\m) = b{t) i+, {o fc }> + Er= h{t) i- {lfc}) ' (3) 

where |{lfe}) is the single-photon state in A:-th mode of 
the reservoir. The amplitude b(t) at any time satisfies 
following integro-differential equation 



b(t) + iu b(t) 



f(t - T)b(r)dT = 0, 



(4) 



where the correlation function is f(t — r) = 
*7 / ^^-e~ iUJk ^ t ~ T ^dk. The non-Markovian memory ef- 
fect of the structured reservoir has been registered self- 
consistcntly in the kernel function f(t— r) in Eq. (4). 
Going back to the double- NVC case, the master equation 
of NVCs can be derived by tracing over the environmen- 
tal degrees of freedom from Eq. (3) [1], 
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(5) 



Hth ^ T^rHOi denote 



where fi(f) = -Im[^] and 7(f) = -Re[^| 
Lamb shifted frequency and decay rate of the NVCs, re- 
spectively. If f{t — t) in Eq. (4) is replaced by a delta 
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function, then Eq. (5) recovers the conventional master 
equation under Born-Markovian approximation [33]. 

Starting from the initial state of the whole system as 
|*(0)) = H-,-)+/3|+,+))|{0 fe }) ri |{0 fe }) r2 [34], we can 
obtain the time-dependent state as 
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.{0k})i\-A0k}h + PW(t))iW)h, (6) 



where Eq. (3) has been rewritten as \ip(t)) = 
&(*) 1+) |0)r + b(t) |-> |l)r with the collective states of 
the reservoir defined as |0) r = |{0fc}) and |l) r = 

^Efc6*(*)|l*> an d b(t) = \/ 1 - \b(t)\ 2 [36]. From Eq. 
(6) the reduced density matrix pNiN 2 f° r the subsystem 
NV1-NV2 can be obtained by tracing over the degrees of 
freedom of the two reservoirs, 



pJViiV 2 (t) = 



f \(3\ 2 \b(t)f (3a*b{t) 2 \ 

p 

Op 

y P*ab*(t) 2 



(7) 
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with p = \pb(t)\ 2 b(t) 2 and q = 1 - \P\ 2 \b{t)\ A - 2p. Simi- 
larly, the corresponding reduced density matrices can be 
obtained for other subsystems like reservoirl-reservoir2 
(riT2) and NVC-reservoir (N\ri, N\T2, N^Ti, A^t^)- 

The quantum correlations can be quantified by EoF 
[29] and QD [28, 37]. The former is defined as E{p) = 

i+yi-c( P ) 
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j with H[x] — —j. i<J52 a/ 
and C(p) = max{0, \/Ai~ — — V^3 
decreasing-order-arranged quantities Ai are the eigenval- 
ues of the matrix p(a£ ® Uy)p*(cfy ® Oy) with p* the 

complex conjugation of p and cr^ B ^ the Pauli matrix 
acting on the subsystem A(B) [38]. On the other hand, 
QD is defined as the minimum difference between two 
ways defining mutual information (MI), Q{p) = I(pab) — 
max{n fe } I(pA\{n.k})- Here I(pab) = S(p A ) + S(p B ) - 
S(pab) is the quantum MI and max{n k } I {PA\{Tlk}) = 
max.[ Uk y[S(pA) - J2kPkS(p A \^k)] is the maximal MI 
when a measurement is performed on subsystem B [28], 
where S{-) is the von Neumann entropy, {11^.} is a com- 
pletely positive-operator-valued measure on the subsys- 
tem B and pu is the respective probabilities. 

Fig. 2 presents the quantum correlation dynamics 
(EoF and QD) of the NV1-NV2 subsystem with differ- 
ent initial states, where two typical cases of ojq < w c and 
cj 3> uj c are considered, corresponding to the NVC's 
transition frequency within and far beyond the band gap 
of the PBG material, respectively. In the former case, 
as shown in Fig. 2(a,b,c), the correlations calculated by 
EoF and QD exhibit obvious oscillation due to the en- 
ergy exchange between the NVC and the memory reser- 
voir, and then the correlations approach a definite value 
in the long-time limit. It implies that the decay rate of 
the excited state \A 2 ) of the NVC approaches zero after 
some oscillations, which leads to the QCP. However, in 
the wo 3* w c case, the residual correlation by EoF and 
QD between NV1 and NV2 vanishes in the long-time 



= —x Ioq 



(1-x) log 2 (l-ar) 
a/A~4}, where the 



limit, although there exists transient revival in EoF [see 
Fig. 2(e,f)]. All the differences between these two cases 
indicate that the PBG plays a vital role in the quantum 
correlation evolution of the NVCs confined in the struc- 
tured environment. We will return to this point later. 
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FIG. 3: (Color online) Time evolution of the quantum cor- 
relations calculated by QD (thick-line) and EoF (thin-line) 
for different subsystems, where the thick-solid line, thick-dot- 
dashed line, thick-dashed line, and thick-dotted line denote 
the QD1, QD2, QD3, and QD4, respectively. The thin-solid 
line, thin-dot-dashed line, thin-dashed line, and thin-dotted 
line denote the Eofl, Eof2, Eof3, and Eof4, respectively, (a) 
a = 1/V2, and ljo = w c /10; (b) a = 0.2, and ljo = uj c /W; (c) 
a = 1/V2, and luo = 10o» c ; (d) a = 0.2, and ujo = 10w c . The 
parameter r\ = 0.2 is used. 

To get a clear picture on how the quantum correlation 
is distributed among different subsystems, we have calcu- 
lated quantum correlations for different bipartite parti- 
tions, such as rcservoirl-reservoir2, NVl-reservoirl, and 
NVl-reservoir2 labeled by QD2 (Eof2), QD3 (EoF3) and 
QD4 (EoF4), respectively, and the quantum correlations 
between NV1-NV2 are denoted by QD1 (EoFl). Fig. 3 
shows the time evolution of EoF and QD in each sub- 
system under different initial states, and one can find 
that for each partition, the correlations denoted by EoF 
and QD behave similarly with obvious oscillations. In- 
terestingly, once the two NVCs are initially prepared in 
maximally entangled state with a = [3 = l/y/2, QD is 
always larger than EoF for both ujq < lj c and loq 3> lo c 
cases, as shown in Fig. 3(a,c). The only exception is the 
curves for QD2 and EoF2 in Fig. 3(c) in the w ^> uj c 
case, which do not obey the monogamic relation. In con- 
trast to this feature, if the NVCs are not initially pre- 
pared in maximally entangled states, as in Fig. 3(b,d), 
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FIG. 4: (Color online) Illustration of the formation of EFBS 
(a) and its dynamical consequence to the excited-state pop- 
ulation (b) when u)a = oj c /W (solid line), ui c (dashed line), 
2ui c (dot-dashed line), 5w c (dot-dot-dashed line), and 10.0o; c 
(dotted line). The intersect of y(E\) — Ei with the horizontal 
axis characterizes the formation of EFBS. r\ = 0.2 is used. 



QD is usually smaller than EoF. Furthermore, it is im- 
portant to note that the QCP in each subsystem exists 
only for the NV's transition frequency within the band 
gap, as shown in Fig. 3(a,b). Another obvious feature is 
that the quantum correlation initially prepared between 
NVCs distributes between other subsystems with differ- 
ent weights. However, for ujq 3> uj c , the residual quan- 
tum correlation of QD1, QD3, and QD4 vanishes in the 
long-time limit, but the quantum correlations (QD2 and 
EoF2) between the two independent reservoirs approach 
their maximum values, which implies that quantum cor- 
relation initially prepared in NVCs has been fully trans- 
ferred to the reservoirs, significantly different from the 
multi-site distribution of quantum correlation in the case 
of Wq < w c . 



IV. THE DISCUSSION AND CONCLUSION 

The Born-Markovian approximation is only valid in 
our case when the correlation function of the radia- 
tion field has an infinitesimal time scale, which is al- 
ways valid when its density of states is a continuous and 
smoothly varying function in the frequency space [39]. 
However, the distribution function for PBG reservoir is 
not a monotonous function due to the existence of band 
edge, where the density of states varies rapidly with a 
rate comparable to the NVCs spontaneous emission rate 
[12], then it becomes non-exponential and the emission 
spectrum becomes non-Lorentzian. By accurately calcu- 
lating the non-Markovian dynamics, wc really find that 
the quantum correlation dynamics displays an oscillatory 
behavior and non-exponential tendency. Figs. 2 and 3 
indicate clearly that the photon emitted through spon- 
taneous emission goes back and forth between the NVCs 
and their local reservoirs. It means that the environ- 
ment is not robust any more during the interaction with 
the quantum system, but changes at the time scale of its 
memory time to reach a new equilibrium state. 



We attribute the mechanism of this novel QCP to 
the interplay between the formed EFBS and the non- 
Markovian effect. Here we argue that there are two es- 
sential conditions for realizing the QCP. One is the exis- 
tence of the EFBS, which provides the ability to preserve 
quantum correlation, and the other is the non-Markovian 
effect, which provides a dynamical way to preserve the 
quantum correlation [40]. Thus, whether the EFBS ex- 
ists or not plays a key role in achieving the QCP, even if 
the reservoir has already shown the memory effect. 

To verify this argument, we find under what condition 
the EFBS can be formed. The total excitation number 
N = (T + cj_ + J2k^k^ k °f Eq. (1) is conserved, so the 
Hilbert space is split into the direct sum of the subspaces 
with definite N. From the eigen equation in the N = 1 
subspace, we have the eigenvalue E\ fulfilling 



-dk = Ei, 



(8) 



where i] = uj 2 d 2 /6ir 2 eqc 3, is a dimcnsionlcss coupling con- 
stant, and 77W0 is the natural spontaneous emission rate 
of the excited state IA2) of NVC. Combining with Eq. 
(2), we can find that y(Ei) is monotonically decreasing 
in the region E\ 6 (— 00, w c ), i.e. in the band gap of the 
reservoir. It means that Eq. (8) may have one and only 
one solution in this region if the system parameters fulfill 
y(uj c ) < uj c [sec Fig. 4(a)]. The curves when wo is smaller 
than or comparable with uj c definitively have an intersec- 
tion point with the horizontal line so that one discrete 
root exists for Eq. (8). We call this discrete eigenstate 
as the EFBS. On the other hand, no root exists in the 
region (w c ,+oo) because it would make y(E\) divergent 
otherwise. As a stationary state of the whole system, the 
EFBS does not lose any quantum coherence during the 
time evolution, implying that the decoherence can be in- 
hibited. Fig. 4(b) shows the dynamical consequence of 
the EFBS on the excited-state population of the NVC. 
We really find that once the EFBS is formed, the excited- 
state population will stabilize to a steady value in the 
long-time limit. 

The underlying physical picture is that during the 
evolution the NVCs excited level \A-i) experiences an 
anomalous "giant" Lamb shift and splits into doublet 
levels only when the NVCs transition frequency u>o lies 
inside or near the band gap [41-43]. One of the doublet 
levels retains in EFBS lying within the gap, whereas the 
other is shifted out of the gap and exhibits resonance flu- 
orescence. In the Wo < w c case, once the FEBS has been 
established, the excited-state population of the NVC in a 
FEBS remains constant in time because the FEBS state 
is actually a stationary state with a vanishing decay rate 
during the time evolution. The claim above can be fur- 
ther verified by the fact that the QCP is always absent 
in the ujq 3> w c case, as shown in the right panels of Fig. 
2 and in the lower panels of Fig. 3, where the EFBS 
does not appear due to the NVCs transition frequency 
far beyond the upper band of the PBG. It is physically 
originated from the fact that the eigenstate of the Hamil- 
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tonian (3) with the eigenvalues resides in a continuous 
band experiences an out-of-phase interference during the 
time evolution, which makes the excited-state population 
approaching zero asymptotically and causes the severe 
decoherence to the NV1-NV2 subsystem. 

In summary, we have investigated the dynamics of 
quantum correlation between two initially entangled 
NVCs embedded, respectively, in PCC. We have shown 
the existence of QCP resulted from both the memory ef- 
fect of the non-Markovian structured reservoir and the 
EFBS. With the mechanism of QCP, we expect to have 
a significant step toward the future full-scale quantum in- 
formation processor based on the increasingly-developed 
nanoscale solid-state technology. It is noted that al- 
though only 2D structure of photonic crystal is taken 
into accounted in our work, our result could potentially 
be applicable to the 3D case with a more complicated 
PBG structure [8]. We expect that the QCP could be 
achieved whenever the EFBS is formed under the condi- 
tion that the transition frequency of the NVC is within 



the band gap of the 3D photonic crystal. The qualitative 
verification to this expectation is worth to be performed 
by combing the explicit form of 3D photonic crystal struc- 
ture in the future. 
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